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We investigate the phase diagram of the compact U{ 1) lattice gauge theory in four dimensions 
using a non-standard action which is invariant under continuous deformations of the plaquette 
angles. Just as for the Wilson action, we find a weakly first order transition, separating a confining 
phase where magnetic monopoles condense, and a Coulomb phase where monopoles are dilute. We 
also find a third phase where monopoles are completely absent. The topological action offers an 
algorithmic advantage for the computation of the free energy. 


I. INTRODUCTION 


In recent years lattice Quantum Field Theory has seen 
a surge of efforts to construct new lattice actions which 
aim at improving the approach to the continuum limit. 
The best-known strategy is that advocated by Symanzik, 
where irrelevant operators of higher and higher dimen¬ 
sion are added to the “standard” (e.g. Wilson plaquette) 
action, with coefficients adjusted perturbatively or non- 
perturbatively to cancel discretization errors of the corre¬ 
sponding power in the lattice spacing a mi- This kind 
of improvement is thus a parametric one, allowing for a 
faster approach to the continuum limit than exhibited by 
the “standard” action. 

However, this is not the only possible strategy for im¬ 
provement. It has long been recognized that departure 
from the continuum limit is more violent for large fields, 
so that suppressing these large fields produces a non- 
parametric improvement [3]- for instance, this happens 
when one trades the Wilson action for the Manton ac¬ 
tion [j|], based on the length of the geodesic in group 
space, or for a “perfect” action [5]: large fields, corre¬ 
sponding to small values of the plaquette trace, are more 
suppressed than with the Wilson action, and at the same 
time continuum behavior is better approximated for a 
given value of the lattice spacing a. 

A more radical suppression of large fields is achieved 
by imposing a strict cutoff: for instance, in a spin model 
one can demand that neighboring spin angles do not dif¬ 
fer by more than a limiting value; or in a gauge theory, 
one may require that the plaquette trace be larger than 
a limiting value. The best-known example of the latter 
is the positive-plaquette action for SU( 2) lattice gauge 
theory [SHE]. While the approach to the continuum limit 
is also improved in this strategy, an important side-effect 
may happen. Localized topological defects can only form 
if the cutoff is not too restrictive. For instance, an 0(2) 
spin model on a square lattice can support vortices only if 
the spins can rotate by 7r/2 or more between neighboring 
sites. If not, the disordered phase of this system disap¬ 
pears entirely. Thus, the cutoff may change the phase di¬ 
agram of the model. A similar situation occurs in lattice 
gauge theory: as pointed out by Liischer (9], if the pla¬ 


quette trace is restricted to “admissible” values greater 
than about 0.97 (for SU( 2)), changes in the topological 
charge become impossible, and topology becomes well de¬ 
fined on the lattice. Topological sectors arise as in the 
continuum theory. 

Here, we consider the extreme strategy where the ac¬ 
tion consists only of a cutoff. In other words, the action 
takes only two values: 0 if all cutoff restrictions are satis¬ 
fied, +00 if not. This kind of action has been called topo¬ 
logical uni, because it does not have any classical small-a 
limit, and the action remains invariant under small ad¬ 
missible deformations of the field. A simple example of 
topological action for an O(N) spin model is: 


S=Y, R o{S i -S j ), Rf) W 

(i,j) 


0 x > cos 6 
Too else 


(1) 


Topological actions raise an interesting puzzle: as the 
constraint between neighboring spins becomes more re¬ 
strictive, the correlation length increases and diverges; 
but what is the action associated with this continuum 
limit? Several studies have investigated different spin 
models m in], and it has been shown in analytically 
solvable O(N) models that the continuum limit is that 
associated with the usual, sigma-model action. In higher 
dimensions numerical investigations also support this 
claim very strongly. 

Here we want to investigate the properties of a topo¬ 
logical action in a gauge theory, and consider the sim¬ 
plest case, namely compact (7(1) lattice gauge theory 
in 4 dimensions. Aside from the continuum limit, we 
also want to study the phase diagram of this system. 
With the Wilson action, a first-order phase transition 
separates a strong-coupling, confining phase and a weak- 
coupling Coulomb phase. This phase transition is as¬ 
sociated with condensation of magnetic monopoles in 
the strong-coupling phase (12]. With a topological ac¬ 
tion, the constraint on the plaquette trace, when restric¬ 
tive enough, is going to make it impossible for magnetic 
monopoles to exist. This may completely alter the phase 
diagram of the theory. 

Finally, topological actions may be interesting for al¬ 
gorithmic reasons: it may be computationally easier to 


2 


move in the space of admissible configurations since they 
all have the same action. While this promise has not yet 
been realized for the Monte Carlo update of such config¬ 
urations, in spin models or in the gauge theory we study, 
we show below that extracting the free energy (or equiva¬ 
lently here, the entropy) is extremely simple numerically, 
and yields valuable information. 

Our paper is organized as follows: we discuss the topo¬ 
logical action of our model in Sec. II, the consequences 
for magnetic monopoles in Sec. Ill, the helicity modulus 
in Sec. IV, propose some arguments about the contin¬ 
uum limit in Sec. V, and discuss how to obtain the free 
energy in Sec. VI. Our results on the phase diagram are 
presented in Sec. VII, followed by conclusions. 

II. THE ACTION 

The obvious analogue of restricting the angles between 
neighboring spins in a spin model is to restrict the real 
part of the trace of each plaquette in a gauge theory. The 
action then depends on one coupling a and is given by 

1 ReTrt/p > a VP 
0 else 

where P denotes a plaquette. Note that this formulation 
is independent of the gauge group but that we from now 
on consider only C7(l) where ReTrUp = cos Op. We could 
thus equally well consider a restriction of the plaquette 
angle Op with \9p mod 27r| < (5 max = arccosa. It is also 
important to note that the link angles, being gauge vari¬ 
ant, are completely unrestricted. The most efficient way 
to generate configurations is to apply heatbath updates 
to the links one at a time under the constraint that no 
plaquette angle exceeds the allowed value. In principle 
this is realized by just uniformly sampling the interval 
[0, 27 t] until an acceptable angle has been found but in 
some cases it might be more efficient to explicitly con¬ 
struct the allowed range of values for the link to be up¬ 
dated. Note that a Metropolis update based on the old 
value may not be ergodic since the admissible region of 
link angles may not be connected. See Fig. [TJ However, 
there are some additional caveats to this kind of single 
link update which will become clear in the discussion of 
the magnetic monopoles. 

III. MAGNETIC MONOPOLES 

An elementary cube on the lattice contains q mag¬ 
netic monopoles if the outward oriented, physical (dp £ 
[—7r,7r[) plaquette angles of its faces sum up to 27t q Ha¬ 
lt is easy to check that q £ {0, ±1, ±2} and that a cube 
with q monopoles must have at least one face with physi¬ 
cal plaquette angle \dp\ > |g|7r/3. This immediately tells 
us that for 5 max < n/3 there cannot be any monopoles 
and the topological action does not describe the same 


^max > tt/2 6 

max < 71-/2 




FIG. 1. Forbidden regions (hatched areas) and allowed re¬ 
gions (black lines) for the angle of a link surrounded by three 
(the others are omitted for clarity) staples Si. When the re¬ 
striction angle S max > n/2 ( left, panel) the region can be dis¬ 
connected whereas if it is smaller than 7t/2 (right panel) it 
will always be connected. <5 m ax is the angle between an arrow 
and the edge of the hatched area opposite to it. 

(lattice) physics as the Wilson action (13]. In fact, a 
change of variables from link to rescaled plaquette angles 
dp/8 max can be used to see that all i5 max < 7r/3 are equiv¬ 
alent up to trivial rescalings. Let us therefore concentrate 
on angles larger than that. 

One might think that if there is a deconfinement tran¬ 
sition at some restriction angle 5 max then it should be 
at <5 m ax = 7 t/ 3 since this angle separates the region of no 
monopoles from a region with monopoles. This turns out 
to be wrong. In a sense this is analogous to the situation 
with the Wilson action. At the deconfinement transi¬ 
tion the monopole density jumps down, but it does not 
jump to zero. The system can sustain a small density of 
monopoles without being confining. The same happens 
for the topological action with a deconfinement transi¬ 
tion at a significantly larger restriction angle than 7r/3. 
Still, there is a non-analyticity in the monopole density 
at 8 max = 7 t/ 3, which we investigate further in Sec. VII 
(see Figs. [9] and [l0| . 

A. Creating monopoles 

To study how the monopoles depend on the cutoff an¬ 
gle <5 m ax it is important to understand what the low¬ 
est monopole excitation is. It is well known that every 
monopole is connected to an anti-monopole via a Dirac 
string and that the monopole worldlines must form closed 
loops on the dual lattice. The shortest such loop has 
four vertices and Euclidean length 2\f2a where a is the 
lattice spacing, and the smallest excitation is thus two 
monopoles and two anti-monopoles each located in one 
of the four cubes sharing a single plaquette. See Fig. [2] 
for an illustration. 

It is also important to consider how such a configura¬ 
tion is created from a configuration with zero monopoles. 
In order to create a monopole in a given cube we need 
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FIG. 2. The smallest possible nontrivial loop of monopoles 
world lines which has Euclidean length 2\/2 a. The y and a 
dimensions are collapsed into one so that each cube is repre¬ 
sented by a plaquette and each plaquette by a link. The fat 
link represents the plaquette shared by all four cubes which 
contain a monopole. A +(—) in a plaquette symbolizes a 
positively(negatively) charged monopole in the corresponding 
cube. 


to change its flux by 2tt at the same time as we respect 
the constraints on the plaquette angles. It is therefore 
relevant to investigate the smallest constraint angle for 
which a change of 27r in the flux is possible. If we update 
a single edge of a cube we will change two of its six pla- 
quettes. The sum of these changes must be 27r and the 
required angles can be minimized by letting the change be 
distributed equally over all involved plaquettes. Hence, 
the restriction on the plaquette angles gives <5 max > 7r/2 
to create a monopole with a single link update. This 
means that for n/2 > J max > 7r/3 the single link update 
is not ergodic and cannot be used on its own. To have an 
ergodic algorithm we need to update at least three faces 
of a cube at the same time, which can only be done by 
updating more than one link at a time, as illustrated in 
Fig-! The minimal update to achieve this is shown in 
the lower part of Fig. [3] where two links of a given pla¬ 
quette are updated together. This update changes three 
plaquettes in each of the four cubes sharing the plaquette 
common to the two updated links, and we thus have a 
chance to create four monopoles down to <5 max = 7r/3 as 
required. 


IV. THE HELICITY MODULUS 

The helicity modulus was first introduced in the 2d 
XY-mo del m where it quantifies the response of the 
system to a twist in the boundary conditions. Because 
the twist is a boundary effect the helicity modulus is an 
order parameter for a system with one massive (finite 
correlation length) and one massless (infinite correlation 
length) phase. This is precisely the case of 4 d lattice U( 1) 
gauge theory where the confining phase features massive 
photons whereas they are massless in the Coulomb phase. 
In the context of a gauge theory the twisted boundary 



FIG. 3. Monopole creation with a single link update (upper 
panel) and a multiple link update (lower panel). The fat 
links are the ones updated and the flux of 7r is spread over 
the plaquettes on the right which means that the single link 
update is ergodic down to 5 ma x = 7r /2 and the two-link one 
to 7 t/3. 


conditions can also be thought of as an external electro¬ 
magnetic flux m- More precisely, we define the helicity 
modulus as 


h = 


d 2 f(4>) 


dcf) 2 


i>=o 


( 3 ) 


where / is the free energy density in the presence of the 
external flux (f>. The flux is introduced by the replace¬ 
ment 


cos (dp) —» cos {Op + <j>) (4) 

for all plaquettes in a given stack of plaquettes, i.e. all 
plaquettes in the set {P^, v (x) \ p, = /io,^ = vq^x^ = 
Xq,x u = x/o} • The orientation and position of the pierced 
stack is arbitrary and with a suitable change of variables 
the flux can also be spread out evenly over the (/zo, U) bi¬ 
planes. For the Wilson action h is a simple difference of 
expectation values 



where the sum in the second term is over all plaquettes 
in the stack defined above. For the topological action on 
the other hand, it is not possible to explicitly perform 
the derivatives. However, since the action for each con¬ 
figuration is the same, the free energy is given solely by 
the entropy, i.e. by the number of configurations with 
a given flux <fi. This can be measured by promoting the 
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flux to a dynamical variable, which is updated along with 
the link angles m By measuring the probability distri¬ 
bution p((f>) (via a histogram method for example) of the 
visited fluxes one thus obtains the full 2n periodic free 
energy m and the helicity modulus 

_ a 2 log p((j>) 

d(j) 2 

Alternatively, and more accurately, one can use all the 
global information from p{(j)) = e~ and fit it to the 
classical ansatz m 

m = - log ^ e -^w - 2 - fc ) 2 

k 

= -logi?3 - ^log27T/3ii, (7) 

where p R plays the role of the renormalized coupling in 
the Coulomb phase and $ 3 (z,q) is a Jacobi theta func¬ 
tion. From this ansatz we can extract the curvature at 
<j> = 0 , i.e. h. analytically and we thus obtain both the he¬ 
licity modulus and the renormalized coupling at the same 
time. We further note that they approach each other ex¬ 
ponentially fast for large ftp. Together with eq. |5]) we see 
that this means that p R m h —> ft— 1/4 as /3 —► 00 , which 
is to say that the coupling constant is not renormalized 
in the continuum limit which is of course common knowl¬ 
edge. 

V. CONTINUUM LIMIT 

It is important to dwell a little on the matter of a 
continuum limit for the topological action. Since all the 
plaquettes are forced to unity when <5 max -A 0 one expects 
that in this limit the correlation length diverges and thus 
that it defines a continuum limit. This point of view was 
examined more thoroughly by Budczies and Zirnbauer 
in m- These authors consider a general weight function 
Wt(Up ), which is a function of a plaquette variable Up 
and some parameter (coupling) t. Granted that there 
exists a t c such that w tc {Up ) = S(U P — id) and that for 
t 7 ^ t c the weight function is some smeared version of the 
J-function, then the lattice gauge theory with partition 
function 

Z t = f d[U] l[w t (U P ) ( 8 ) 

J P 

has a continuum limit as t —> t c . Furthermore, the au¬ 
thors claim that under “ favorable conditions” , the con¬ 
tinuum theory will be Yang-Mills theory. It is not pre¬ 
cisely defined what conditions are considered favorable, 
but close to the identity element, the plaquette variable 
is well approximated by Up = e la ~ Fp ~ 1 +ia 2 Fp—a i Fp. 
Thus, in order for the continuum action to be oc f Tr F 2 
the weight function w t certainly has to satisfy some con¬ 
ditions on the moments of the tangent vectors of the Lie 


group. At the very least the first moment must vanish 
and the second moment needs to exist and have the cor¬ 
rect sign. The authors indeed give an example in [16] of 
a weight function, in two dimensions and for gauge group 
U(iV), which satisfies the <5-function condition but which 
has the wrong continuum limit. The problem is identi¬ 
fied with the non-existence of the second moment for the 
considered weight function. 

The topological action which we use clearly satisfies the 
(5-function constraint since the weight function has sup¬ 
port only on a compact region of width ±(5 max around the 
identity element and thus goes to S(Up — id) as (5 max —> 0. 
Because of the compact support and invariance under 
Hermitian conjugation we also conclude that the first 
moment vanishes and that the second is positive as it 
should. It is therefore probable that this action will have 
the correct quantum continuum limit and indeed all nu¬ 
merical evidence suggests that it does. 

A simple check one can perform is to use for wt(Up) 
a combination of angle restriction and Wilson plaquette 
term with negative /?. By taking (5 max —► 0 the action still 
satisfies the (5-function constraint but the negative value 
of P will try to bend the distribution in the wrong direc¬ 
tion to make the second moment of Wt negative. Clearly, 
for a fixed value of (j the action will still be almost flat 
as long as (5 max is small enough, so in order to change the 
continuum limit, ft needs to be taken to —00 at the same 
time as <5 max —» 0. Then, if the magnitude of /3 is large 
enough we expect that the continuum limit is spoiled. 
This can also be observed in numerical simulations, and 
although it is somewhat of a pathological example it still 
gives some insight as to when one can expect to obtain 
the correct continuum limit. 


VI. FREE ENERGY 

Here, we show how to evaluate the free energy, analyti¬ 
cally in a Id toy model, and numerically for more realistic 
cases. 


A. Id XY model 

Consider a periodic chain of N spins s,; £ 0(2) with 
a topological action which restricts the angle of each 
link ti = Sisj +1 to be smaller than (5 max . Let U = 
exp (iOi), 9i £ [— <5 max ,<5 max ]. The partition function of 
this model then takes a very simple form, 

(5max tv / / N \ \ 

(9) 

and describes a collection of N non-interacting, con¬ 
strained links with the only condition that the product of 
all links is one. The normalization of the angle integrals 
serves to keep Z finite as the number of links is taken 
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to infinity and is just a subtraction of the ground state 
energy. 

The total angle can take values 27rm, m £ 
{— [ jV 2^ ax J )■••> [ — 2^ x J } and thus to is the winding 
number or topological charge of the system. The par¬ 
tition function can be expressed solely in terms of the 
total angle by convoluting the uniform distributions of 
the individual links N times. The distribution of the 
sum of N i.i.d. uniform variables converges very rapidly 
to the normal distribution, in this case with zero mean 
and variance N6^ nax / 3. Anticipating the N —> oo limit 
we thus neglect the small deviations from the normal dis¬ 
tribution and write 


Z = 


2nNSl. 


NS„ 


' dd exp - 


3(9 2 


2 NS 2 . 


S (exp ( id ) — 1 ) 


-ns b 


V3 

(2n ) 3 / 2 a 


[y/N q;J 

exp 

m=— [\/iVaJ 



( 10 ) 


where we have defined a = v/ ^ < ^ I * ax . We can now take 
N —> oo whilst keeping a fixed to obtain 


the two representations by 

exp (—(exp (— ) 

Xt_ ^a(ex P (-^)) 

1 / _ exp (— 27 ) ^3 ( ex p(~^)) \ 

■^l 1 #3 (exp (-^)) ) 


(14) 


Since the elliptic function and its derivative are analytic 
functions V/3 £ R + there is no phase transition but there 
are two distinct regimes with a rather abrupt crossover. 
In the low temperature regime, (3/1 > 10, the partition 
function is almost independent of (3 and the topological 
susceptibility is very close to its zero temperature value 
(47 t 2 /) -1 whereas in the high temperature region, (i/I < 
10, the partition function is approximately \/2n(3 and Xt 
rapidly drops to zero. 

Note that, when NS max < 2i r, topological excitations 
are forbidden and \t = 0. However, the continuum limit 
is obtained while keeping (NS^^ fixed, so that the lat¬ 
tice spacing varies oc J 2 lax . Therefore, in this Id model 
the parameter region where Xt = 0 disappears in the 
continuum limit. 


Z = 


V3 

(27r) 3 / 2 a 


OO 

51 exp 

m =—oo 



(27r) 3 / 2 a 


t?3 




( 11 ) 


where $ 3 ( 3 ) = $ 3 ( 0 , q) is the third Jacobi elliptic theta 
function. Since the sum in the partition function is over 
the winding number to it is straightforward to calculate 
(to 2 ) and the topological susceptibility Xt = j (to 2 ). In 
the limit /3 = Na ot a 2 oo (where a is the lattice spac¬ 
ing) one should find Xt = where I is the moment of 
inertia of the quantum rotor which the model describes. 
This allows us to determine a in terms of /3 and I and 


the result is a = 



which leads to 


Z = 




( 12 ) 


With Poisson’s summation formula we can go from the 
winding number representation to the energy representa¬ 
tion in which 


Z = ti 3 





(13) 


It is now evident that the excited states are doubly de- 

>2 

generate and the energy differences are E k — E 0 = §y as 
is well known. The topological susceptibility is given in 


B. Higher dimensions and gauge theories 

In higher dimensions, due to the lattice Bianchi iden¬ 
tities, the integration over the constrained variables no 
longer factorizes and we can not calculate the partition 
function analytically anymore. However, in the small 
<5max regime where there are no topological defects the 
partition function must be 

Z = (2J max )" dof (15) 

(or one, depending on the normalization), where rid.o.f 
is the number of independent degrees of freedom. As 
the topological defects are turned on, the functional de¬ 
pendence on <5 max will change and there will be a high 
order and practically undetectable phase transition. As 
<5max is further increased the topological defects will start 
to play a more important role and eventually the real 
phase transition of the model will occur. If one would 
have access to the partition function, or free energy, one 
could directly extract the properties of the transition. 
Fortunately, since the topological action is constant, the 
partition function is pure entropy and can thus be mea¬ 
sured by Monte Carlo simulations by simply counting 
the number of configurations at a given value of <5 ma x- 
If Figure [2] we show the derivative of the free energy 
density / = — V~ 1 logZ with respect to <5 max for the 
2d XF-model (upper panel) and the 4 d U(l) gauge the¬ 
ory (lower panel) for various lattice volumes, obtained by 
Monte Carlo simulations. It is clear that the derivative 
is smooth in the XT-model where the transition is of 
infinite order (BKT) and that it is discontinuous in the 
1 /( 1 ) case where the transition is first order. 
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COS <f m ax 



FIG. 4. The derivative df /d5 m ax of the free energy density 
/ = — V -1 log Z for the 2d XT-model (upper panel) and the 
4 d U( 1) gauge theory (lower panel ) as obtained from Monte 
Carlo simulations with a topological action. There is a clear 
distinction between the smooth derivative in the XT-model 
which has an infinite order phase transition and the discontin¬ 
uous behavior, signaling a first order transition, in the U( 1) 
gauge theory. The vertical line marks the critical restriction 
cos 5 C . In the case of the XY-model it has been taken from HD 
where it was extracted from a fit of the diverging correlation 
length. 


VII. RESULTS 


Let us now turn to the numerical results. Primarily 
what we are interested in is the phase structure of the 
model and the order of the possible deconfinement transi¬ 
tion. To this end we have measured the monopole density 
and the helicity modulus as a function of the restriction 
cos6 max . We compare these results with the correspond¬ 
ing observables obtained with the Wilson action in Figs. [5] 
and [6] it is obvious that the transition is even weaker 
than the weak first order transition seen with the Wilson 
action. We can try to quantify the strength of the transi¬ 
tion by fitting the helicity modulus in the confining phase 


using a simple model of a first order transition PHE]. 

h= l + X- 1 exp(-UA f(x-x c )Y (16) 

where h + is the helicity modulus in the Coulomb phase 
(which is assumed to be constant), A/ is the latent heat, 
X is an anisotropy factor between the two phases and 
x is the coupling, either (3 or cos<5 max . After taking fi¬ 
nite size effects into account the best fit is shown as the 
lines in Fig. [6] The data is well described by the ansatz 
and one finds that the fitted value of the latent heat for 
the Topological action is about half of what it is using 
the Wilson action, which is consistent with the weaker 
transition seen in the monopole density. 

To further establish that the transition really is first 
order we show histograms of the monopole density close 
to the transition for three different volumes in Fig. [7] A 
double peak structure is formed and enhanced as the vol¬ 
ume increases, which is a clear indication that the tran¬ 
sition is first order. Also the Monte Carlo history shows 
clear tunneling events between two metastable states. 
Together with the discontinuity in the first derivative of 
the free energy with respect to the cutoff we conclude 
that the topological action has a first order transition at 
^max ~ 1.95. 

To determine the characteristics of the two phases 
we look at how Wilson loops of different sizes behave. 
Naively, we expect an area law when i5 max is close to 
7 T since the interaction between plaquettes will be very 
weak, as for the Wilson action where j3 < 1. This can 
be seen in the left panel of Fig. [lj If the forbidden re¬ 
gions become very narrow then the individual links are 
hardly influenced by their neighbors and each plaquette 
angle is more or less uniformly distributed in the interval 
[—<5 max , <5max] which gives an average plaquette trace of 
sin(<5 max ) /A max . For a loop with area A, this is raised to 
the A’th power. For restrictions <5 max close to zero on 
the other hand, the links are heavily influenced by their 
neighbors (right panel of Fig. [T]) and the total angle of 
the loop should depend on the perimeter rather than the 
area. This is demonstrated in Fig. [8] where we show the 
Creutz ratios 


x(R) 


(W(R,R))(W(R-1,R-1)) 
(W(R,R-1)){W(R-1,R))’ 1 ’ 


where W(I, J) is a planar, rectangular Wilson loop with 
sides I and J. We have performed the R —» oo extrapola¬ 
tion under the assumption that the corrections are of the 
form e~ R . Note that this is not a precise measurement 
of the string tension but rather a characterization of the 
two phases. We have also checked that the magnitude of 
the Polyakov loop acquires a vacuum expectation value 
in the low monopole density phase. 

Another interesting thing to investigate is how the 
monopole density depends on the renormalized coupling. 
The monopole mass is proportional to (3r = e]\ 2 and the 
density decreases exponentially with the mass. This is 
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FIG. 5. The monopole density n for the topological action 
(upper panel) and the Wilson action (lower panel). For the 
Wilson action the first order nature of the transition is rather 
evident even for a 12 4 lattice whereas for the topological ac¬ 
tion we have to go to much larger lattices to see a fairly dis¬ 
tinct jump. 


FIG. 6. The helicity modulus h for the topological action (up¬ 
per panel) and the Wilson action (lower panel). The lines are 
the best fit to eq. (161, which describes the data in the con¬ 
fining phase (the model assumes a constant h in the Coulomb 
phase) very well for both actions. 


a statement about physics so it gives us a direct way to 
compare the two actions. In Fig. [9] we show the monopole 
density as a function of the renormalized coupling and we 
see a clear exponential decay as expected. For the topo¬ 
logical action the decay is significantly faster, which could 
be interpreted as a reduction in the discretization errors: 
for a given effective coupling, there are fewer lattice ar¬ 
tifacts (monopoles) that disturb the order of the system. 
For (j max < 7t/3 the density is even strictly zero and the 
model is completely insensitive (up to trivial rescalings) 
to further reduction of 5 max - 

With a mix of single- and two-link updates we have 
been able to measure the monopole density down to 
densities around 10~ 8 . The exponential dependence on 
cos(i max persists to ^ max ~ 1.69 after which the density 
smoothly changes into a power law in ( 1/2 — cos<5 max ) 
with an exponent which is fitted to be 11.70(6) as can 
be seen in Fig. [To] We tentatively ascribe this change 


of functional behavior to the approach of a phase tran¬ 
sition. A naive argument, which works well in the 2d 
XY -model, leads to a monopole density which is polyno¬ 
mial in the small deviation (i5 max — 7 r/3). The argument 
is based on convolutions of (near) uniform plaquette or 
link distributions. To create a single vortex in the spin 
model close to the threshold n/2 we need to convolve 
the link angle distribution four times, which makes the 
joint distribution oc (4<5 max — 9) 3 for the cumulative an¬ 
gle 9 around a plaquette. This needs to be evaluated 
at 9 = 2n (one vortex) which gives a vortex probability 
oc (<5 m ax — 7 t/2 ) 3 . Vortices always come in pairs so we 
expect that the density is proportional to (<5 max — 7t/2) 6 
which is in good agreement with what we have obtained 
from Monte Carlo simulations. By a similar argument 
one would expect a monopole density oc (i5 max — 7 r/ 3) 20 
due to six plaquettes in 4 cubes containing a monopole. 
The deviation in the power law from the predicted 20 


























Rescaled monopole density 


COS ffuxx 


0.038 

0.036 

0.034 

0.032 
e 

0.03 
0.028 
0.026 
0.024 

0 20000 40000 60000 80000 100000 

Monte Carlo time 

FIG. 7. The probability distribution of the monopole density 
close to the transition for various volumes (upper panel) and 
the corresponding Monte Carlo history for the 24 4 volume 
(lower panel). The distributions are rescaled in such a way 
that one peak is at -1 and the other at +1. That way the 
separation of the peaks in relation to the widths can be di¬ 
rectly compared between different volumes. It is evident that 
the two peaks become more distinct for larger volumes which 
indicates a first order transition. Also the obvious tunneling 
between two different states in the Monte Carlo history backs 
up this statement. 



FIG. 8. The Creutz ratios y(.R) given by eq. (171 as well 
as an R —> oo extrapolation assuming corrections of the form 
e~ R as a function of the restriction cos <5 max obtained on an 8 4 
lattice. There is a clear transition between a confining phase 
with nonzero string tension and a deconhned phase with an 
perimeter law for the Wilson loops. 
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FIG. 9. The monopole density in the Coulomb phase as a 
function of the renormalized coupling fin, for the topological 
and the Wilson actions. The different rates of decay could be 
attributed to different lattice spacings for the two actions. 
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to the observed ss 12 is rather large, but the argument 
does not take into account that the 4 monopoles are not 
independent of each other, so it is not so surprising that 
one finds a smaller exponent. 


VIII. CONCLUSIONS 

We have simulated U( 1) lattice gauge theory using an 
unconventional “topological” action. We find that this 
action describe the same physics as the Wilson action, 
i.e. there is a confining strong coupling phase where 
magnetic monopoles condense and Wilson loops follow 


an area law, separated by a (weak) first order transition 
from a Coulomb phase with an exponentially suppressed 
monopole density and a perimeter law for the Wilson 
loops. We have, in this specific case, not found any con¬ 
crete advantages which would motivate the choice of this 
action over the Wilson action although at a given value 
of the effective coupling in the Coulomb phase there are 
significantly fewer monopoles (lattice artifacts). This is 
in line with other known cases where a topological ac¬ 
tion reduces discretization errors HU. Perhaps the most 
interesting approach is to search for optimized combina¬ 
tions of a standard action and constrained fields. An 
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FIG. 10. The monopole density for the topological action as 
a function of the restriction. There seems to be a smooth 
transition from an exponential decay to a power law at 
COS <5max 0.12. 


interesting feature of the topological action is the direct 
access to the free energy itself. 

One interesting open question is the nature of the extra 
transition at <5 max = 7 t/ 3 where there is a non-analyticity 
in the monopole density as it goes from nonzero to strictly 
zero. A similar phenomenon occurs at <5 max = 7r/2 for 
an XY model, and when plaquettes become restricted 
by the “admissibility condition” in gauge theories. One 
may argue, however, in the (7(1) case at least, that this 
transition will have no impact on the physics because 
the monopole density close to the transition is extremely 
small anyway. 
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